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ABSTRACT 

We measure the velocity dispersions of six, galactic globular clusters using spatially integrated 
spectra, to test for the effects of internal dynamical evolution in the stellar mass-to-light ratios, T*, 
of star clusters. In particular, we revisit whether the low values of T* we found in our previous study, 
from which we concluded that there are at least two population of stellar clusters with distinct stellar 
initial mass functions, are artificially depressed by relaxation driven mass loss. The combination of our 
previous sample of five old clusters and these six now provides an order of magnitude range in cluster 
mass with which to explore this issue. We find no relationship between cluster mass, or relaxation 
time, and Y*. Because relaxation is mass dependent, we conclude that the values of T* for these 
clusters are not strongly affected by dynamical effects, and so confirm the presence of the population 
of clusters with low T*. 

Subject headings: stars: formation, luminosity function, mass function; galaxies: fundamental param- 
eters, evolution 
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1. INTRODUCTION 

A variety of independent lines of evidence now point 
to a stellar initial mass function (IMF) that can vary 
from one environment to the next. In elliptical galax- 
ies the study of po pulation-depende nt spectral fea tures 
(|van Dokkum fc Conrovll201fl; ISpiniello et al.ll2012T> and 
detailed dynamical modeling (|Cappellari et alTl2012D 
suggest varying IMFs that are diff erent than those mea - 
sured in the local Galactic disk ([Bastian et al.1 l2010h . 
in the Small Magellanic Cloud deep star counts sug- 
gest an IMF that does not turn over at low masses 
(|Kalirai et al.l |2012|) , and in the study of stellar clus- 
ters anomalies are found both among M31 clusters 
(jStrader et al.l 1201 lft and in those of the Milky Way 
and its satellites (jZaritskv et al.ll2012l hereafter Paper 
I). While it is necessary to continue to explore the IMF 
in all of these environments and provide more evidence 
in support of these initial findings, the study of clusters 
presents one significant advantage over that in the other 
environments. Because star clusters are the most likely 
to contain only a simple stellar population, of a single 
age and metallicity, these systems provide the possibility 
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of a relatively direct approach to unraveling the cause of 
IMF variations — if the variations are confirmed. 

In Paper I we found a bifurcation in the apparent prop- 
erties of star clusters. Specifically, we found that younger 
clusters (age < 10 Gyr) typically, but not exclusively, 
have discretely larger values of the stellar mass-to-light 
ratio, T», for their age than do the older clusters. The 
explanation we proposed is that there exist (at least) two 
populations of clusters that have different IMFs. There 
are, of course, a number of alternate explanation for the 
difference in properties, which we discussed in Paper I, 
but here we present new data and an empirical explo- 
ration of one of the principal potential sources of system- 
atic error, the impact of internal dynamical evolution on 
T 

Older clusters, which have suffered the greatest degree 
of two-body relaxation, may have preferentially lost large 
numbers of low-mass stars, and thereby have lower T* 
than models of stellar evolution would predict for a given 
IMF. As such, one wonders if the discordant properties 
of the younger and older clusters can be reconciled. In 
Paper I we explored this possibility quantitatively using 
published models o f stellar cluster dynamical evolution 
(jAnders et al.ll2009l ) , but the application of those models 
has large uncertainties because certain input parameters 
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are not well constrained. Here we present observations 
of six additional old clusters, which more than doubles 
our sample of older clusters observed and analyzed in 
the exact same manner. Furthermore, these new clusters 
now enable us to probe the effects of dynamical evolution 
over a larger range of cluster mass. Because the effects of 
two-body relaxation are most pronounced in lower mass 
clusters, a sample with a large mass range provides an 
empirical test of the effects of two-body relaxation. We 
present the data in §2, discuss the derived masses, T* 
and examine the dependence of T* on mass and relax- 
ation time in §3. We find no dependence, demonstrating 
that two-body relaxation has had a minimal impact on 
our results and that the values derived for T* for our old 
clusters are not artificially deflated. As such, we con- 
firm the existence of the lower T* population of clusters. 
We have therefore addressed one of the most significant 
source of systematic uncertainty in our previous results. 
Future work will present more data related to the popu- 
lation of stellar clusters with large values of T*. 

2. THE DATA 

Our spectroscopic data come from a set of obser- 
vations taken with the Las Campanas du Pont tele- 
scope (100-inch), p rimarily for the pur p ose o f chemical 
abundance studies (IColucci et al.l 120111 l2012f) . We use 
the compilation of iMcLaughlin fc van der Marell (2005) 
for the necessary ancillary data (age, half light ra- 
dius, luminosity, and metallicity) . When model data 
are used from the compilation, we choose res ults ob- 
tained using the Wilson models dWilsonl H975I ). which 
McLaughl in fc van der Mare] (|2005| ) demonstrate are su- 
perior in fitting the radial surface brightness profiles of 
these clusters. All of these choices follow our procedure 
in Paper I and, as then, our conclusions are insensitive 
to the choice of the Wilson model fits. 

We obtained our spectra with the echelle spectrograph 
at the 100-inch at Las Campanas (du Pont) during dark 
time in 2000 and 2001. To obtain integrated spectra of 
the clusters, we utilize t he same spect r oscopi c drift scan 
technique described by iColucci et al.l (|2011l) . To sum- 
marize, we set the telescope in motion to raster scan 
the slit across the cluster during the exposure, defining 
both the angular length and height of the raster (both 
set to the same number). The exposure time then sets 
the rate of the scan, such that the full scan is com- 
pleted within the allotted exposure time. The echelle 
slit is l"x 4", allowing uniform coverage of a 32 x 32 
arcsec 2 high-surface brightness region of the cluster. We 
took multiple exposures to homogenize the scanned clus- 
ter region and for cosmic ray removal. The spectra 
have a wavelength coverage of approximately 3700 to 
7800 A, with declining sensitivity and spectral resolu- 
tion toward the blue end. We reduced the spectra us- 
ing standard IRAF routines (see Colucci et al. 2011), 
includi ng the scattered light sub t ractio n described in de- 
tail in [McWiniimSBirnitei^ dHH) . Further details 
of the obs ervations and data reduction can be found in 
IColucci et all (poll . 

The spectra are of the same high quality as in our orig- 
inal study (Paper I) because the total exposure times 
are also typically around 10,000 seconds (see Tabled]). 
Examples of similar spectra and the spectral fitting de- 
scribed below are given in the original paper. There 



are no substantive differences between those data and 
that described here. This similarity is confirmed by the 
internal velocity uncertainties derived from the fitting, 
described below, that are quantitatively comparable to 
those obtained for the clusters in Paper I. The S/N of 
the spectra varies among the objects and as a function 
of wavelength but generally exceeds 20. 

2.1. Measuring Velocity Dispersions 

A discussion motivating our approach is presented in 
Paper I, and we briefly review it here. 

We fit a Gaussian broadening function to each of the 
same set of lines selected in Paper I. We reject any lines in 
the object spectra that are either clearly blends or suffer 
some other complication and, after fitting, we reject lines 
that do not produce an acceptable fit, where acceptable is 
defined by xt < 2.3. To calculate xt we adopt a per pixel 
uncertainty determined from the fluctuations about a flat 
continuum in line-free areas of the spectrum. However, 
we adopt the same uncertainty value for the full spectral 
range for any given cluster. Our results are not highly 
sensitive to this value because we only use these xt values 
to remove questionable fits from further consideration, 
and visual inspection confirms that those lines that have 
been rejected by this criteria are poorly fit. 

Once the line fitting is complete, we calculate the over- 
all best value of a using the same approach as in Paper 
I, where we place different weights on downward and up- 
ward uncertainties because larger apparent broadenings 
can be caused by various systematic issues (blended lines, 
poor continuum fit, focus errors). The final value of a 
for each cluster is calculated using the average of all the 
measurements and x 2 /iV, where N is the number of data 
points, with uncertainties in any individual measurement 
defined as resulting in x 2 /A = 1. To down weight mea- 
surements that are inflated by blends or focus errors, 
we set upward uncertainties to be 3 times larger than 
downward ones. The uncertainty on our final "mean" a 
is derived by identifying the range of c's that generate 
xt~Xnu min < 2.71, corresponding to the 90% confidence 
interval (see Table [lj . For a more complete description 
and discussion of the issues involved we refer the reader 
to Paper I. 

For the subset of our clusters that have been observed 
previously, we compare our a measurements to pre- 
existing values. W e obtain previous measuremen t s from 
the compilation of McLaug hlin fc van der Marell (|2005l) , 
which includes literature values of a for 9 of our clus- 
ters (over the e ntire age range), and other sources that 
include 3 more jlMevlan et al.lll99lt iFischer et~aT1 119931 : 
lLane et al.lfe OlOD. For N GC 362, a val ue of the disper- 
sion was not presented by IFischer et al.l ([T993) so we cal- 
culate the standard deviation of the individual stellar 
velocities they provide, after excluding the sources they 
identified as either non-members or binary stars. The ex- 
ternal comparison is positive in that 9 of 12 (75%) of the 
measurements agree to la and 11 of 12 (92%) agree to 
2cr (Figure [T]). As we discussed in Paper I, the published 
dispersion of the one significant outlier, NGC 1866, is 
highl y sensitive to the inclusion or exclusion of a single 
star (Lane et alJl2010f >. 

3. DETERMINING MASSES 
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Fig. 1 . — A comparison of our measurement of the cluster velocity 
dispersion, a, to previously published measurements by indepen- 
dent investigators. The comparison includes clusters of all ages 
(from Paper I and this study). The line represents the 1:1 rela- 
tion. The agreement is excellent, but our new values provide the 
improved precision necessary to uncover subtle differences in T*. 

In Paper I we applied the p rescription presented 
and tested by iWalker et al.l (|2009f ) to derive masses for 
spheroidal galaxies. As discussed previously in more de- 
tail, our finding that their prescription, when expressed 
in the formalism of galaxy scaling relations, 



log r h = 2 log a - log I h - log T h - 0.73, 



(1) 



where rh, Ih, and T\ are the half light radius, the sur- 
face brightness within that radius, and the mass-to-light 
ration within that radius, almost exactly matches the 
empirical finding of a scaling relation that applies to all 
stellar systems from stellar clusters to massive ellipticals 

log r h = 2 log g - log I h - log T h - 0.75. (2) 

(|Zaritskv et al.l 120081 120 111) provides additional confi- 
dence in this method. The empirical results verify that 
there is little scatter (~ 0.1 dex) about this relationship 
for objects ranging from globular clusters to massive el- 
liptical galaxies. We also demonstrated in Paper I that 
this method results in values of T\ in good agreement 
with published values calculated using more detailed dy- 
namical models. Using Equation^ we evaluate Th for 
the six clusters of this study and present the results in 
Table [T] For systems without dark matter, which we pre- 



sume includes these clusters, T* = T/j. All photometric 
quantities are presented for the V band. The uncertain- 
ties in Y* are calculated using only the uncertainty in a, 
as the internal uncertainties on the other parameters are 
proportionally much smaller. 

3.1. Dynamical Evolution of Clusters and Its Impact on 

T 

A range of dynamical process affect stellar clusters. 
Of particular relevance to the type of study we are con- 
ducting is the evaporation of low-mass stars via two- 
body interactions, which will depress the value of T, and 
has the strongest im pact on the oldest clusters (|Spitzerl 
1958tlK ruiisscn 2008). In Paper I, we used the models of 
Anders et al.l (j2009f ) who calculated the disruption time 
for clusters, ta, which depend on the cluster mass and 
local (external) tidal field. Applying the models, how- 
ever, depends on the choice of the normalization factor £4, 
which is the age at which a 10 4 M Q cluster loses 95% of 
its mass under the corresponding conditions. Although 
each cluster should be assigned its own value of £4, cor- 
responding on its internal density, orbit, and tidal field 
strength, we had to adopt general values of £4 without 
much guidance. Nevertheless, we proceeded to calculate 
the correction factors for Y*, but found that in general 
these were small and even in the most extreme cases 
they were less than a factor of two, which is insufficient 
to reconcile the differences between the two populations 
of clusters. However, because of the uncertainty in the 
selection of £4 and the dependence on idealized models, 
an unsatisfactory level of uncertainty remains. 

The data described here, in combination with that in 
Paper I, provides an opportunity to reexamine the ques- 
tion in a more empirical manner. The modeled relaxation 
effects are inversely proportional to the mass of the clus- 
ter, so we expect low mass clusters to have suff ered the 
most. Specifically in the lAnders et al.l (|2009l ) models 
td oc Af 062 . With the new data presented here, the sam- 
ple of old clusters we have has more than doubled and 
spans an order of magnitude in mass, suggesting disrup- 
tion times that should differ by a factor of 4. Therefore, 
if some of our oldest clusters, age ~ 13 Gyr, have lived 
for a significant fraction of their disruption time, then 
the most massive in the sample have disruption times 
that are > 50 Gyr, and so should show negligible effects 
of relaxation. We therefore expect that if cluster val- 
ues of T* are affected by relaxation, the effect should be 
significantly larger at the low mass end of our range. 

In Figure [2 we plot the relationship between T* and 
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Fig. 2. — T t vs. cluster mass and relaxation time. Relaxation 
time is presented in arbitrary units. We find no discernible trend 
between T* and cluster mass or relaxation time. The lack of trends 
argues against the importance of two-body relaxation effects be- 
cause in such a scenario T* should be lower in the low mass and/or 
short relaxation time systems. Clusters from Paper I are plotted 
as black squares, while the clusters presented here are plotted as 
red circles. 

cluster mass. We find no evident mass dependence. This 
result does not rule out dynamical evolution in all sys- 
tems, for example the two clusters with lower than av- 
erage T* may have experienced some. The variations in 
Y* seen among the older clusters in other studies, and 
presumably i n ours, can be expl ained by two-body relax- 
ation effects (|Leigh et al.l I2012I ). and our results do not 
contradict this claim. Our conclusion is that dynamical 
evolution cannot result in the much larger effect required 
to reconcile these observed values of T* with those of our 
other cluster population (Paper I). The only possible way 
to avoid this conclusion is to evoke a detailed inverse cor- 
respondence between the strength of tidal effects and the 
mass of the clusters (with tidal effects being stronger on 
the the more massive clusters and canceling out the mass 
dependence of the two-body relaxation effects) . We see 
such a correspondence as contrived, particularly because 
some of the older clusters (from Paper I) are in the Large 
Magellanic Cloud, which should provide a significantly 
different tidal environment than for those clusters in the 
Milky Way. 

In Figure [U we also plot T* against a more direct 
estimate of the relaxation time, t r ex Nt cross , where N 
is the number of stars and is proportional to M* and 
tcross oc r^/a. Again, we find no dependence that leads 
us to conclude that our previous estimates of T* were 
biased low by dynamical evolution. 

Finally, in Figure [3] we reprise Figure 10 from Paper I. 
The plot shows the value of T* that a cluster will have 
(or had) at an age of 10 Gyr, T* j io. We describe this cal- 
culation in Paper I, although for the clusters presented 
here this correction is minor because they are all only 
slightly older than 10 Gyr. We have made several al- 
terations from the version of the plot shown in Paper I. 
First, we include now only clusters for which we have 
measured the velocity dispersions (previously we had 
augmented our set with velocity dispersions presented in 



iMcLaughlin fc van der Marel (2005)). This choice pro- 
vides a degree of homogeneity to the data that facilitate 
comparisons among clusters. Second, we have removed 
the correction for dynamical evolution because we now 
consider it to be minor, for these clusters, and highly 
uncertain. Third, we recalculate the uncertainty regions 
for the lower T^io population including these new data 
and excluding the previously included literature data. 
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Fig. 3. — Evidence for two initial mass functions revisited. We 
plot T*io, the value of T* at 10 Gyr for each cluster, obtained 
as described in Paper I, versus iron abundance and age. In red we 
highlight the clusters presented in this paper. The shaded regions 
represent the standard deviations from the mean for that corre- 
sponding population. These are different than shown in Paper I 
because they have been recalculated from the data presented in 
this Figure. 



The distinction between the two populations continues 
to be clear. As noted in Paper I based on literature data, 
but now confirmed with our data, there is overlap among 
the populations along the metallicity axis. Metallicity 
alone, therefore, cannot be the cause for different IMFs. 
Additional data that we are now collecting and process- 
ing will attempt to expand the sample at intermediate 
ages, which will help address the question of whether 
there is an analogous overlap of populations along the 
age axis. 

We close now by evaluating the gains obtained by 
our re-observations of clusters that already have velocity 
dispersion measurements in the literature. The origi- 
nal reason for observing those clusters was to establish 
that there is no systematic difference in the velocity dis- 
persions obtained through our integrated light measure- 
ments and those in the literature based on individual 
stars. Figure Q] confirms that there is no difference. Fur- 
ther examination of that Figure also highlights that we 
are claiming significantly improved precision in our mea- 
surement of the velocity dispersions. This improvement, 
if it can be confirmed, is critical because T has a a 1 
dependence and because among the literature measure- 
ments of the required parameters a has the largest frac- 
tional uncertainty. As such, the uncertainty in a domi- 
nates the uncertainty in T. Without independent mea- 



surements of a that are of comparable or better preci- 
sion, we cannot directly confirm our uncertainty esti- 
mates. However, indirectly we find that the scatter in 
the calculated T^io values for the ancient clusters ob- 
tained with literature ct's (Figure 10, Paper I) is nearly 5 
times larger than that for our sample (Figure |3|), suggest- 
ing that the dispersion in our a measurements is indeed 
significantly smaller that that of previous measurements 
and that further observations of clusters along the lines 
presented here, even if they have existing literature val- 
ues of <7, is of value. 



4. SUMMARY 

We have measured the velocity dispersions of six, 
galactic globular clusters using spatially integrated spec- 
tra, to test for the effects of internal dynamical evolution 
in the stellar mass-to- light ratios, T„, of star clusters. 
We conclude, based on the lack of any detectable varia- 
tion in Y* with cluster mass, that dynamical relaxation 



is not affecting the entire population of clusters. This 
finding addresses one of the principal potential sources of 
systematic uncertainty in our previous argument for dis- 
tinctly different IMFs among two sets of stellar clusters 
in the Milky Way galaxy and its satellites (Zaritsky ct al. 
2012). Because additional concerns remain, and the 
possibility of a non-universal IMF is sufficiently impor- 
tant, we are obtaining more data on other clusters to 
further test our original claims. In the meantime, evi- 
dence for variations of the IMF is appearing in a num- 
ber of environments that i s coming from a wide range of 
independent o bserv ations (Ivan Dokkum fc Conrovll2010l 
Strader et aTll2011t ICaippellari et alj|2012t iKalirai et al 



20121) . We advocate increased caution when interpreting 



extragalactic observations using a single IMF. 
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